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Abstract 



^ ' Recently Pawlowski and Raczka (P&R) proposed a unified model for the fundamental 

interactions which does not contain a physical Higgs field. The gravitational field equation of 
their model is rederived under heavy use of the computer algebra system Mathematica and 
its package MathTensor. 

o . 

^ ■ 1 Introduction 

The computer algebra system Mathematica, together with its package MathTensor, is very useful 



for executing differential geometric calculations on a four-dimensional Riemannian spacetime, as 



was shown by Soleng for example. Here we want to demonstrate the effectiveness of these tools 
by picking an example from the current literature which is of some fundamental importance. 

The gravitational sector of a unified model for the fundamental interactions, proposed by P&R 
| H , is determined by the Lagrangian density 

£geom = V^9 (1 + P) - A($t$) 2 - pC ijk l , (1) 

with $ as Higgs field and R (curvature scalar) and Cyw (conformal Weyl curvature tensor) as 
gravitational fields depending on the metric tensor gij of the Riemannian spacetime, g := det(<7y). 
The coupling constants /3, A, and p are dimensionless because of the conformal invariance of the 
model. 

The total Lagrangian density of P&R, 

^gcom ~\~ Afield ? (2) 



*To appear in Hehl, F.W., Puntigam, R.A., Ruder, H. (eds.): Relativity and Scientific Computing - Computer 
Algebra, Numerics, Visualization. Lectures given at the 152 th WE-Heraeus Seminar, 18-22 Sept 1995, Bad Honnef 
(Germany). Springer, Berlin, 1996 
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depends additionally on the field part, which subsumes the contributions from all nongeometrical 
pieces. The explicit form of £fi c id is of no relevance to usfj] since it is absorbed in the definition of 
the 'material' energy-momentum tensor 

T ij ._ _2_ Afield , 3 s 
Hamilton's principle yields the gravitational field equation 

^gcom ^-^field 

8gij 5gij 

or 



= , (4) 



1 <5£gcom _ 1 rpij ^ 



V~9 s 9ij 

The total Lagrangian (^) is conformally invariant. Therefore it is possible, within the P&R model, 
to fix a conformal gauge for the scalar Higgs-type field $ according to 

v 2 

<F $ = — = constant . (6) 



Our goal will be the computation of the left-hand side of (g) by means of the computer algebra 
tools mentioned above. Thereby we want to check the corresponding results of P&R 0. 



2 Riemann Tensor and its Irreducible Pieces 

Before we can commence with our calculations, we have to get hold of the definition of the Weyl 
curvature tensor Cijki- It is known that the curvature tensor of a four-dimensional Riemannian 
space has three irreducible pieces: the Weyl tensor Cijki, the tracefree Ricci tensor fCij, and the 
curvature scalar R. 

We take the conventions of Misner et al. || and define the curvature tensor {V % jk = Christoffel 
symbols) 

R l jki '■= dtT l ji — diT l jk + r J m fe r m j-j — r%„/ v m jk , (7) 

the Ricci tensor and the curvature scalar 

Rij ■= R k tkj and R := g ij R tj , (8) 
respectively, and, eventually, the tracefree Ricci tensor 

jtif= Rij - - g^ R. (9) 

The metric has signature ( — h ++)• 

The Riemann tensor (20 independent components) decomposes into three irreducible pieces, 
20= 10© 9© 1, or 

Rijki = {l ) Rijki + (2) R llk i + (3) R llk i , (10) 

jFrom the Higgs-type Lagrangian, two pieces feature in the geometrical Lagrangian, namely the first two terms 
in (hi), whereas the rest of it is attributed here to the field Lagrangian in (H). 
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with the definitions 

( 2 )% w := >H, I fc 8 l l and ^R tJ kl := i *J . (11) 

The first irreducible piece is traceless and has to be identified with the Weyl piece: C«jy := ^'Rijki ■ 
If resolved with respect to Cijfc/, Eq. ([lO]) can be read as the defining equation for the Weyl tensor. 

Let us use MathTensor in order to perform the irreducible decomposition defined in ( |l0| ) and 
(pT|). MathTensor recognizes TraceFreeRicciR. Therefore the pieces ^Rjjki, denoted by Rl, R2, 
and R3, can be computed by an almost verbatim translation of ( JTc| ) and ([IT]) into MathTensor: 



Dimension=4; Rcsign=l; (*Default*) 

Def Unique [TraceFreeRicciR [la_ , lb_] , ,PairQ [la, lb] ] 

R2 [li_ , lj_ , lk_ , 11_] =2 Antisymmetrize [Antisymmetrize [ 

TraceFreeRicciR [lk, li] Metricg[lj ,11] ,{li,lj>] ,{lk,ll>] 
R3 [li_ , lj_ , lk_ , 11_] =Expand [Antisymmetrize [ 

1/6 ScalarR Metricg [lk, li] Metricg [lj , 11] , {li , lj>] ] 
Rl [li_ , 1 j _ , lk_ , 11_] =RiemannR [li , 1 j , lk , 11] 

- R2[li,lj ,1k, 11] - R3[li,lj ,1k, 11] 



We shall verify this decomposition with MathTensor. We first use the predefined tensor WeylC in 
order to check if the definitions are correct. 

Furthermore we want to make sure that our decomposition is really irreducible. For this purpose 
we compute the traces of ^Rijki , and find - as well as the usual symmetries of the curvature tensor 
- that Cijki is traceless, the Ricci tensor of ^Rijki is JRfy, and ^Rijki has neither a traceless piece 
nor a traceless Ricci piece, rather only the curvature scalar with the correct factor 1. The input 
reads: 



TraceFreeRicciR [li , lj] / . TraceFreeRicciToRicciRule 

dif f l=CanAll [ (WeylC [li , 1 j , lk , 11] -Rl [li , 1 j , lk , 11] ) 

/ . WeylToRiemannRule] 
dif f 2=CanAll [Tsimplif y [CanAll [dif f 1 

/ . TraceFreeRicciToRicciRule] ] ] 

TraceR2=Expand [Metricg [ua , ub] R2 [la , li , lb , 1 j ] ] 
TraceR3=Expand [R3 [la, li , lb, lj] Metricg [ua,ub] ] 
TraceRl=Tsimplif y [Expand [Metricg [ua,ub] Rl [la, li , lb, lj] 
/ . TraceFreeRicciToRicciRule] ] 
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3 The Topological Euler Density 

Let us come back to (0). Its last term is proportional to 

C 2 :=sf=-gC ijkl C iikl . (12) 
Densities will be denoted by script letters. The computation of the Bach tensor [|| 

1 SC 2 



(see also [|5|) can be simplified, if one splits off a divergence term from (|12|). Such a term is the 
topological Euler density || 



£ := V^9£ abcd R ll abR kl cd e ljkl . (14) 



We can show that ([14|) represents a divergence: 

f = a *y^ 4 ) (is) 

with 

V 1 := V=9£ ijkl £ab cd T a CJ (]-R b dkl + \ T b mk T m d ) . (16) 



2 3 

Furthermore, also by means of MathTensor, the explicit form of ( p"4j ) turns out to be 

8 ;= IS ( Rl3kl Rm ~ 4 RlJ ^ +R2) - (17) 
The corresponding input reads: 



Eulerdl= 

-1/(128 Pi~2) Epsilon[ua,ub,uc,ud] Epsilon [li , lj ,1k, 11] 
RiemannR[ui ,uj ,1a, lb] RiemannR [uk,ul , lc , Id] 
Eulerd2=CanAll [Expand [Eulerdl/ . EpsilonProductTensorRule] ] 

We substitute into (|l2] ) the definition of the Weyl tensor and eliminate the emerging curvature 
square piece by means of (|l7|). Then we find the simplified formula 

C 2 = C 2 + d l V l , (18) 

with 

C 2 := 2 ( /,"• /,',, - ^) = 2 {0* Mij ^) , (19) 
which has also been cross-checked by means of MathTensor: 
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WeylSquarel=WeylC[li,lj ,1k, 11] WeylC[ui,uj ,uk,ul] 
WeylSquare2=CanAll [Dum [WeylSquarel/ .WeylToRiemannRule] ] 
checkCsquare=Expand[WeylSquare2 - 32 Pi~2 Eulerd2] 
checkCsquare2=Expand [checkCsquare 

/ .RicciToTraceFreeRicciRule] 



4 Bach Tensor 

In order to compute the variation of the remaining piece C 2 of the Lagrangian, we follow closely the 
scheme that was demonstrated in the lecture by Soleng Q. After initalization, we use Variation 
in order to compute 5C 2 , followed by a series of partial integrations, rules, and simplifications. 
The last step actually computes the variational derivative: 



varCl=Variation [Sqrt [-Detg] *checkCsquare ,Metricg] 
varC2=PIntegrate [varCl .Metricg] 
varC3=PIntegrate [varC2 ,Metricg] 

varC4=Canonicalize [Absorbg [ApplyRules [varC3 , RiemannRules] ] ] 
Bachl [ui_ ,uj_] =Tsimplif y [VariationalDerivative [Expand [ 
varC4/Sqrt [-Detg] ] , Metricg , li , 1 j] ] 

MetricgFlag=True 



The outcome of this variation reads (a semicolon denotes the covariant derivative) 

2 1 
Bachl = Bij = - /,':,:, - 2 R tr y k + - 9lj R. k ;fe 

— —R 2 gij + — RR^ + R k i R kl — AR kl Rikji , (20) 

B^ = Bji and = B k k = . (21) 

The result © differs slightly from that of P&R.f] 

In order to compute the Bach tensor (|20|), we could have used (|l^) rather than (|l9|). In this 
case ( |20| ) would pick up two additional terms that compensate each other, as is explicitly shown 
in We feel, however, that the present detour, via C 2 , pays off in conceptual and computational 
simplicity. 

2 Their result reads, see Q Eq. (7.2)] and Eq. (5.2)]: 

2 -k 2 -k 

~ g + 2 Rij;k' ~ — 9ij R;k' 

~3 R 2 9rj + 2 RR H + 9n Rki R kl -4R kl Rikji ■ 

The third term of the first line carries an incorrect factor two. Therefore their Bach tensor is no longer traceless, as 
is required by conformal invariance. Up to a (conventional?) sign, their first line is identical to that of (|2o[), whereas 
the second lines coincide. 
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5 The Bach Tensor Streamlined 



The Bach tensor takes on a more transparent form if we express the curvature pieces in ( |2C| ) 
exclusively in terms of the irreducible pieces. 

By default, MathTcnsor docs not recognize that the trace of C^u vanishes over abritrary 
indices. Therefore we explicitly define rules to take care of this fact. After this preparatory step, 
we can directly reformulate the Bach tensor. Thus we put in: 



Def Unique [WeylC[la_,lb_,lc_,ld_] ,0, PairQ [la, lb] I I 

PairQ[la,lc] II PairQ [la, Id] II PairQ[lb,lc] II 
PairQ[lb,ld] II PairQ [lc , Id] ] 

Bach2 [li_ , 1 j _] =Tsimplif y [CanAll [Expand [Bachl [li , 1 j ] / . 

RiemannToWeylRule/ .RicciToTraceFreeRicciRule] ] ] 



This computation yields 

2 It ■ k 

Bach2 = Bij = — R.^-j — g-Rjfc' 9ij ~ 2 f^ij;k' 

+-^Rflij +^J$ki Mj —gijflkiM —^ft Cikji . (22) 

If we introduce 

7^=^?"^^, (23) 

where 

7fo .«= 7 yM=0 and f r>-,/' (). (24) 
then the tracelessness B k k = and the symmetry B^ = Bji of the Bach tensor become manifest: 

2 

+ ^ RJtij +4 7 y kl Mmk M"\ - 4M H C ikjl . (25) 

The trick for the corresponding computation is again the use of a series of rules that expresses the 
Bach tensor in terms of kl (denoted by Gam) : 



Def ineTensor [Gam, {{2 , 1 , 3 , 4} , 1 , {1 , 2 , 4 , 3} , 1}] 
Def Unique [Gam [li_ , 1 j _ , lk_ , 11_] , , PairQ [li , 1 j ] ] 

RuleUnique [GamRulel ,Metricg [li_ , lj_] 

TraceFreeRicciR[lk_,ll_] TraceFreeRicciR[lm_,ln_] , 
(- 4 Gam [li , lj ,um,uo] 

+ 4 Symmetrize [Metricg [li , urn] Metricg [lj ,uo] , {um,uo}] ) 
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TraceFreeRicciR[lm,ln] TraceFreeRicciR [lo ,un] , 
PairQ [lm , lk] &&PairQ [11 , In] ] 
RuleUnique [GamRule2,Metricg[li_,lj_] CD[ScalarR,lk_,ll_] , 
(-4 Gam[li,lj ,ul,uf]+4 Symmetrize [Metricg [li ,ul] 
Metricg[lj ,uf] ,{ul,uf>] ) CD [ScalarR, 11 , If] , PairQ [lk, 11] ] 
Bach.3 [li_ , 1 j _] =Tsimplif y [Dum [Expand [Bach2 [li , 1 j ] 

/ . GamRulel/ . GamRule2] ] ] 
TraceBach3= Expand [Metricg [ui ,uj] Bach3 [li , lj] ] 



6 Gravitational Field Equation of the P&R Model 

The gravitational field equation (||) can now be made explicit by substituting (Q), (^J), and ( |l5| ) 
into it: 



(l + f3)v 2 1 S(y^R) \v 4 1 5y/=g~ 
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'9 &9i 



^9 Sgi 



' 2 



We compute (by MathTensor) the variations (G u = Einstein tensor) 

1 6{y/=g~R) 1 



^9 % 



G lJ 



Rv + -Rg v 
2 



and 



1 5^—g 1 y 

= t: Q 

2 y 



(26) 



(27) 



(28) 



V~9 °9ij 

By inserting these relations into (|2^), we find the (corrrected) P&R field equation (see J?], (Eq. 7.2)] 
and §, (Eq. 5.2)]) 

(l + (3)v 2 „„-, Xv 4 1 
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■G y + g l] +pB 11 



Provided /3 ^ — 1 and v 2 0, we can put (E9) in a more conventional form, 



Gi 



3Xv 2 



' y 2(l+/?) 9ij (l+/3>2^- (i + /3)t,2"«' 
or, after some (computer) algebra (T := T fe fc , ?"ij' = Ty — T gij/A): 



Yip 
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■ li. 



3Au 2 12p 6/1 

y ' + 2(1 + /?) " (1 + /3)«2 By _ ~(l + /3) v 2 ^ _ 2 Tffu 



(29) 



(30) 



(31) 



We can decompose this equation into its two irreducible pieces, the tracefree and the trace piece: 

v 2 v ~ fl + /3V^ y ' 

(32) 



'•' (1 + /3K y 
6Aw 2 



(1 + /3> 2 
6 



(1 + /3) 



(1 + /3)w 2 



T. 



These two pieces ( |32|) of the (corrected) gravitational field equation ( |2^ ) or (||l]) of P&R, together 
with the explicit form of the Bach tensor ((2^), represent the general result of our considerations. 
Note that in the trace piece of ( |32| ) there occur only second derivatives of the metric, in contrast to 
the fourth order derivatives featuring in the Bach tensor of the tracefree piece. Incidentally, models 
leading to somewhat similar gravitational field equations have been discussed since the early 1920s 
by many people (see || |f| and pi , and the literature cited therein). 



In vacuo we have 



Yip 



( 33 ) 



„ _ 6\v _ . 

u ~ ~TTp ~ : 4Acosm - 

A glance at (|25| ) shows that we can find a special solution of the vacuum equation ( |33| ) by using 
the Einstein vacuum equation with a suitable cosmological constant as an ansatz: 

Rij = A cosm Qij (34) 

or, alternatively, 

I *l - "iA < 35 > 



Then the Bach tensor ( |25| ) vanishes, By = 0, and (|33j) is fulfilled. 

We have collected the MathTensor code that verifies the results of this section: 



Def ineTensor [B , {{1 , 2} , 1}] 
Def ineTensor [T,{{1 ,2} , 1}] 
Def ineTensor [TrFrT , {{ 1 , 2} , 1}] 
Def Unique [B [li_ , 1 j _] , , PairQ [li , 1 j ] ] 
Def Unique [T [li_,lj_] ,T, PairQ [li , 1 j ] ] 
lagHE=Sqrt [-Detg] ScalarR 
varHEl=Variation [lagHE,Metricg] 
varHE2=PIntegrate [varHEl ,Metricg] 

HE [ui_ ,uj_] =Expand [1/Sqrt [-Detg] VariationalDerivative [ 

varHE2,Metricg,li,lj]] 
lagconst=Sqrt [-Detg] 

var const l=Variation [lagconst , Metricg] 

const [ui_ ,uj_] =Expand[l/Sqrt [-Detg] VariationalDerivative [ 

var const 1 , Metricg , li , 1 j] ] 
MetricgFlag=True 

FieldEql [li_,lj_]=(l/12 (1+beta) v"2 HE[li,lj]) + 
(1/4 lambda v~4 const [li , lj] ) + 
(rho B[li,lj]) - (1/2 T[li,lj]) 

FieldEq2 [li_ , 1 j _] =Simplif y [Collect [Expand [-12/ ( (1+beta) v~2) 
*FieldEql [li,lj]] ,{RicciR[li,lj] , ScalarR}]] 
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TraceFieldEq2=Simplif y [Expand [Metricg [ui ,uj] 

FieldEq2 [li.lj]]] 

RuleUnique [TRule,T[li_,lj_] , 

TrFrT[li,lj]+l/4 T Metricg [li.lj]] 
FieldEq3 [li_ , lj_] =Tsimplif y [Expand [FieldEq2 [li , lj] - 
1/2 TraceFieldEq2 Metricg [li , lj] ] ] 
TraceFieldEq3=Expand [Metricg [ui ,uj] FieldEq3 [li , lj] ] 
TraceFreeFE [li_ , 1 j _] =Expand [FieldEq3 [li , 1 j ] - 

1/4 Metricg [li,lj] TraceFieldEq3 
/ .RicciToTraceFreeRicciRule/ . TRule] 



7 Discussion 

Models very similar in their gravitational sectors to that by P&R have been developed, amongst 
others, by Gregorash & Papini |ll], [l2| and in ]l3|, 14| , the latter one, though, in a metric-affinc 



spacetime with additional conformal invariancc (lj, Sect. 6]. Since conformal invariancc is ac- 
commodated much more naturally in a spacetime with a Weyl piece, we believe that these post- 
Riemannian model should be reconsidered in the light of the more recent developments. 
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